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We investigate the Compton scattering vertex of charged scalars and photons in scalar quantum 
electrodynamics (SQED). We carry out its non perturbative construction consistent with Ward- 
Fradkin-Green-Takahashi identity (WFGTI) which relates 3-point vertices to the 4-point ones. There 
is an undetermined part which is transverse to one or both the external photons, and needs to be 
evaluated through perturbation theory. We present in detail how the transverse part at the 1-loop 
order can be evaluated for completely general kinematics of momenta involved in covariant gauges 
and dimensions. This involves the calculation of genuine 4-point functions with three massive prop- 
agators, the most non-trivial integrals reported in this paper. We also discuss possible applications 
of our results. 

PACS numbers: ll.15.Tk, 12.20-m, ll.30.Rd 

I. INTRODUCTION 

Evaluating Green functions for gauge theories in the non perturbative regime is as important as it is difficult. For 
example, understanding confinement and dynamical mass generation in quantum chromodynamics (QCD) requires 
knowledge of quark and gluon propagators in the infra red. These objects are tied to the 3- and 4-point vertices through 
Schwinger-Dyson equations (SDEs) as well as Slavnov- Taylor identities (STI), For example, unlike abelian gauge 
theories, 4-point ghost-ghost-quark-quark scattering kernel arises even at the level where these identities relate 2- 
and 3-point vertices in covariant gauges. 4-point quark-quark scattering kernel for off-shell external legs is another 
important quantity. Its infra red behaviour in the coordinate space would tell us how the strong potential rises 
for large distances between two static quarks. Study of the 4-gluon vertex can improve our insight of the running 
coupling in the infra red region, Understanding non perturbative behaviour of a 4-point function is a daunting 
task, so a detailed study of such a function in a relatively simpler scalar quantum electrodynamics (SQED) in terms of 
Ward-Fradkin-Green-Takahashi identity (WFGTI) Q can provide an important first step towards much more involved 
similar functions in QCD. 

We start by providing the most general basis in terms of which we can expand out the complete Compton scattering 
vertex. We use an appropriate set of basis vectors so that the components of the 4-point vertex longitudinal to the 
external photons can be readily identified. By invoking WFGTI, which relates the 4-point Compton scattering vertex 
to the 3-point scalar-photon vertex and the 2-point scalar propagator, we construct the longitudinal vertex to all 
orders in perturbation theory in terms of lower n-point functions. This WGTI-conserving construction will enable us 
to go to the next order of approximation in SDE studies. In principle, we shall now be able to truncate it at the level 
of the 4-point functions rather than the 3-point functions. Note that the WFGTI leaves the transverse part of the 
vertex undetermined which has to be evaluated using the brute force of perturbation theory. 

Despite the fact that perturbative evaluation of n-point Green functions in gauge theories is receiving more attention, 
analytical results for arbitrary gauge and dimensions with completely off-shell external legs even for 1-loop 3-point 
vertices for SQED and QCD P,H, H, 0] have been reported not so long ago. Results for spinor QED can be derived 
by appropriate replacement of color factors in the quark-gluon vertex. These perturbative results provide a natural 
guide to the possible non perturbative structures of the fermion-boson vertex, structures which are vital in the reliable 
truncations of Schwinger-Dyson equations at the level of the 3-point functions to study confinement, dynamical mass 
generation and hadron spectrum, see for example (sl. [ol. [l0| . 

A natural next step in this direction is the evaluation of a 4-point function in arbitrary gauge and dimensions 
for off-shell external legs. To our knowledge, no such complete calculation exists for any gauge theory though a 
lot of work exists which deals with the calculation of box-diagrams in particular kinematical regimes related to the 
problem at hand. For example, a method to calculate massive Feynman integrals using Mellin-Barnes technique was 
developed in pT| for propagators and triangles with arbitrary powers of the denominators. The first on shell infra red 
divergent scalar box was calculated in [l^l using a mass regulator for an internal photon; this result has extensively 
been employed in many electroweak and heavy quark calculations since then. Photon-photon scattering box diagrams 
were treated in [isj using different representations for 1-loop integrals. 2-, 3- and 4-point massive integrals for 
arbitrary momenta and dimensions were studied in with explicit examples of the integrals relevant to Bhabha 
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scattering. The 4-point master integral in D dimensions using McUin-Barnes technique was studied in llSll . General 
techniques for solving n-point scalar massive and massless 1-loop integrals were developed in [H, In |l8l|. the box 
diagram in arbitrary dimension was calculated for on-shell external legs with internal propagators of varying masses 
in connection with the first order radiative corrections to Bhabha scattering in D dimensions. 1-loop amplitudes for 
4-point functions with two external massive quarks and two external massless partons were studied in detail in p^ . 
In connection with 2 — > 2 QCD scattering amplitudes, scalar 1-loop integrals in different kincmatical regimes have 
been evaluated in [13, [U, [22|. More recent work on the subject can be found in [2^. Our present article extends 
these earlier endeavours to all external legs being off-shell with completely different degrees of off-shellncss. Moreover, 
instead of calculating one particular topology, we evaluate a complete process in a gauge theory (Compton scattering 
in SQED) in an arbitrary covariant gauge, calculating all the diagrams contributing to the process in general number 
of dimensions D. 

SQED is a simple gauge theory as compared to spinor gauge theories (in the sense that there is no Dirac matrix 
structure), but the loop integrals have the same basic topologies as the ones that appear in spinor QED and QCD. This 
makes it an attractive theory. Recently, analytical result for the six-photon helicity amplitudes has been reported for 
SQED, [23|. Apart from the fact that we develop the machinery needed to obtain the off-shell 1-loop 4-point function, 
the study of processes such as Compton scattering 

e-(p)+7(fc)^e-(p')+7(fc') (1) 

has a value in its own right. For example, the electric and magnetic polarizabilities of a pion are related to the 
Compton scattering amplitude at threshold. Due to the fact that the internal structure of the pion contributes little 
to these polarizabilities, 1-loop results are important to calculate (see [2^ and references therein) because these play a 
key role in the calibration of high energy colliders. Moreover, in terms of physical insight, the analysis of processes with 
virtual or real photons attached to an electron or a quark line help us probe the electron or nucleon substructure (see 
[2^ and references therein for a review in such developments). In particular, the calculation of 1-loop amplitudes with 
external massive quarks are crucial to describe heavy quark hadroproduction and constant theoretical developments 
have led to numerical or semi-analytical procedures that provide results with improved precision. These advances have 
allowed for a continuous flow of different mathematical methods enabling scalar and tensor loop integral evaluation 
with massive lines . 

In this article, in addition to constructing a WFGTI conserving longitudinal Compton scattering vertex, we also 
present its complete 1-loop calculation for off-shell external momenta and for arbitrary gauge and dimensions for 
massive SQED. This involves 2-, 3- and 4-point scalar and tensor integrals up to 2 indices ^. The box integrals 
involved have at most three massive propagators. We proceed by writing down the contributing diagrams and 
calculating them in terms of scalar and tensor integrals which we list in the appendices. Due to exchange symmetries, 
only eleven of the twenty eight diagrams are independent. All the calculations have been carried out using FORM 
or/and Mathematica 6.0. 

We have organized the article as follows: In section II, we start out by proposing a convenient basis to decompose the 
Compton scattering vertex into its longitudinal and transverse parts. Moreover, we construct the longitudinal vertex 
non perturbatively. It exhausts 5 of the 10 basis tensors. In Section II, after setting the notation and introducing 
the Feynman rules, we present all the 1-loop topologies for the Compton scattering and write them in terms of tensor 
and scalar integrals to be evaluated. The reduction of the tensor integrals in terms of scalar ones and their evaluation 
has been done in the appendices. Subtraction of the longitudinal part earlier constructed yields the undetermined 
transverse vertex to the one loop order. This should provide us with a guide towards its possible non perturbative 
extensions. In section IV, we present our conclusions. The appendices have been dedicated to the calculation of the 
necessary integrals and presenting the 1-loop results for the full 4-point vertex. 

II. WGTI CONSERVING LONGITUDINAL VERTEX 

The WFGTI which relates 3- and 4-point Green functions is : 

k'^T,^{p',k'-p,k) ^ T,{p + k,p)-T,{p',p' ~k) , 

k''T,^{p',k'-p,k) ^ T,{p',p' + k')~T,{p^k',p) , (2) 
where the vertices involved are full, as shown in the diagram below : 



Tensor integrals involving three Lorcntz indices also arise but the two diagrams involving them cancel each other out because of their 
charge conjugation symmetry. 
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The first argument of the 3-point vertex throughout represents the momentum of the incoming fermion and the second 
that of the outgoing one. The WFGTI identities allow us to write "longitudinal parts" of the 4-point vertex in term of 
the 3-point vertex. These parts are defined as the ones which vanish on carrying out any of the contractions described 
in Eq. ^ . The remaining "transverse parts" will have to be calculated by the brute force of perturbation theory. We 
conveniently define 

= k,{p + p')-k' -k-k'{p + p% , 
Rf.1 — • k ~ k k^^ , 

f^jy — ^ ' k k kjj . ( ^ } 

Note that this definition is suitable to separate out the transverse part of the vertex from the longitudinal one. Thus 
the complete 4-point vertex can be written in its most general form as : 



rj. = Ag^. + Bn{k ■ k'g,^ - k,k'^) + BriQ^k'^^ + B^^R'.k'^ 
+ B2ik^Q^ + B22Q i,Q i_i + B2zRyQ 

+ B^ikuRfj, + B32Q'^Rp. + BssR'j^R^ , 



(4) 



where the coefficients A, B12, B21, -B31 and B13 are determined by WFGT identity to all orders in perturbation 
theory. Therefore, 



r^i/ — A9a'1' + BuQ'jyk' + BisR^k' 



- fj,. 



B2 



Bn{k ■ k'g,^ - k,k') + B22Q',Q,. + B23<Qm + B32Q',R^ + Bs^R^Rf, . 



(5) 
(6) 



The coefficients which constitute the longitudinal part of the vertex are 
1 



A 
B21 
B12 
B31 
B13 



k ■ k 



- {A-i(p + k) - A-^ip) + A~Hp' ^k)- A-i(p')} 



1 rA-i(p + fc) - A-i(p) A-i(p') - A-i(y-fc) 



k- k' [ {p + kY - p2 p/2 _ (^pi _ ^•)2 

1 / A-i (y - fc') - A-^y ) A-i (p) - A-i (p + fc') 
k-k' 
1 



{k-k'Y 
1 



{p' — k'Y — p'^ p^ — [p + k\ 

{[{p + kf - p']Tt{p + k,p) + [{p' - kf - p'']Vt{p',p' ~ k)} 



-k^[rT{p + k,p)-TT{p',p' ~k)] 

-k'^iTrip' -k',p')-rT{p,p + k')] 



{[{p' - kr-p"]rT{p' - k\p') + [{p + kr-P^]rT{p,P + k')} 



[k-k'Y 

where Tt is the transverse part of the 3-point scalar-photon vertex defined as 

A-i(p + fc)-A-i(p) 



T^{p + k,p) - (2p + fc), 



{p + kY — P^ 



2(fcpP -k-k Pf^)TT{p + k,p) . 



(7) 



(8) 



Note that in all those calculations, where only the longitudinal projection of the true vertex is involved, Eqs. ([7]) 
constitute the exact non perturbative result. However, Bn, B22, B23, -B32 and i?33 remain undetermined and need 
to be calculated perturbatively order by order. 

In the next section, we calculate the transverse vertex to one loop order in perturbation theory. This requires the 
calculation of the full vertex to that order. 
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III. COMPTON SCATTERING VERTEX 
A. Preliminaries 

In this sub-section, we shall set the notation and define all the scalar and tensor integrals to be evaluated. We 
define the bare quantities in the usual form: the scalar propagator A(p) = — m?), the photon propagator Aj^^, ^ 

~ [gfj,^p^ ^ (1 ~ OPfj-Pi^] Iv^i the 3-point vertex r° — (k and the 4-point double photon vertex e^r"^, — e^gf^^, 

where ^ is the general covariant gauge parameter (such that ^ = corresponds to Landau gauge) and e is the usual 
QED coupling constant. Following is the list of 1-, 2-, 3- and 4-point scalar integrals that we need to evaluate : 

T = / T-^—T. - (9) 



^i^) = / ^lul V2 21 ' (10) 

^ ' ^ J K -to2][(/c-u;)2 -m2] ' (^^) 

= / 4Ul '^''Z 21 ' (12) 

f d^w 

J w'^[{k — w)'^ — m'^][{p — w)^ — m^] ' ^ ^ 

^^^'^^ " / K - m2] [{k - w)2 - m2] [{p - w)2 - m2] ^^^^ 

J(fc,p) = / , (15) 

' J ^^[(fc — i/;)2 — m2][(p — ui)2 — m2] ' 

^ '^''^) J w'^[{k — w)'^ — in'^][{p — w)'^ — m'^][{q ~ w)-^ — m'^] ' ^ 



We also meet tensor integrals up to two indices in such a way that the numerator is or w^w'^ with the same 
denominators as in Eqs. ([51 (17)) . Therefore, we shall employ the same notation as above for these tensor integrals 
with the only difference that we shall add an appropriate tensor superscript (or subscript) to the representing symbol. 
Thus, for example, 

I^'^ik,p) = f ^"^'^ d^'^ . (18) 

This set of integrals appear in all the possible 1-loop Feynman graphs contributing to the Compton scattering process, 
as we shall see next. Results of all these integrals are given in detail in the appendices. 

Some of these integrals diverge in 4-dimensions and need to be regularized and renormalizcd. We do not carry out 
renormalization in this article but it is worth mentioning that the scalar integrals T, K{k), L{k), J{k,p) and V{k,p) 
are divergent, whereas, I{k,p), I{k,p) and U{k,p,q) are convergent. As the tensor integrals can be decomposed in 
terms of scalar integrals, it is easy to deduce their convergence properties. 
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In this section, we express all these diagrams in terms of the scalar and multi-indexed integrals defined in the 
previous section. Note that we use the notation q = k + p and q' = k' + p throughout : 



^d,Ap, q,p') 



^ D20 



^'i^,-{k,p,q,p') 



2ie/ 
2 



■ 9 
(2^ 



■{-g^''K{q) + {l-OL^Aq)} , 



(g) + ^\{ K^q) + (1 - Oip^ + q^L^q) 



{q^ — m?) 



[q^ — 711?) 



g^'^l Kip - p') - 2{p + p'ri^{p,p') + 4p ■ p'Iip,p') 



(2^) 



+ (e - l)[K{p-p') - 2{p + p'ri>.{p,p') + ^p'^P^'J^uip.p')] 
]^{2{p + qTp^'Iip, q)-ip + qrrip, q) - ^p'l^ip, q) + 2/^'^(p, q) 

-(1 - [{p^ - + qT^ip. q) - 2(p' - 9) 

+2L^%q)~{p + qrL''iq)]} , 



ie^ (p + qY 



^D^Ap^q^p') = 



-ie^ [p' + q'Y{p + q'Y 



(p + qYip' + qr 



{ - 2p"^K{p') + K'^ip') + (1 - 0{P'^ - m^)LAp')} 
2 ^nY \ , T 



(Y — rn? ' [q^ — m?Y 



K{q)- 



(1-0 



L{q) 



2{2tt)D 
K{q) 



[cY — ir?Y q^ ~ 

- Aq'^K{q') -T + Aq'^K^{q') + (1 - e) [T - Aq'^K'^iY) 
+Aq'JpL-Aq')\^, 

K{Y-q) , K{p') 



^ 2 - 1 Mb , 9) - , 



+ 



Y — m? 



-2ip + qY 



9 9 ? 

AY -q , - p'^ .\Ti^( I ^ 
q'^ — rw^ (7^ — 777^ 



9 9 9 
m"^ q'^ — 



Y" + q" 
2(g2 — m? 



■K'iY - q) 



K'^iY) , K-{q) 



9 9 9 9 



/9 9 

— 



-(1-0 



{p + qYiY + qY 



[Y^ -m^).J{Y,q)-L{q) 



9 9 



i(p') 



-2(p + g)^ 



q^ — m'^ 



-le 



^{ (fc + 2p)'' (fc + p + p')" [(-2m2 + (p' - p)2 - 2p • p')C/(P, + /(? - P, y - P) 



(19) 
(20) 



(21) 



(22) 
(23) 
(24) 



(25) 



(26) 



2(27r; 

-I{p, q) - /(p', q)] -2{k+p + YY{{-2Tr? + (p' - pY - 2p ■ Y)U^{p,Y . q) 

+P^{q -p,Y-Y- IAp. q) - IAp'. q)] - (2fc + 4p)^[(-2m2 + (p' - pY - 2p ■ Y)U''{p,Y, q) 

+Jy -p,Y-p)- i'Yp, q) - iAp, q)] + m-^m^ + ip' - p? - 2p • p')u^Ap,p', q) 

+/'"'(g - p,p' - p) - /^"^(p, q) - /'"-(p', q)] + [{k + 2p)^(fc + p + YY [(-™'(p' + Y^ - m") 

+P^Y^)v{p,Y, q) + (p' - 2m^)u{p,Y, q) + i{q - p.p' -Y + (™' - p')J{p. q) 

+ {iiY - p")J{Y, q) - I{P, q) - I{p\ q) + Hq)] - 2(fc + p + p'Y [(-™'(p' + Y' ~ riY) 
+p^Y^)V>'ip,Y,q) + (p' +p" - 2m^)U^'{p,Y,q)+I''{q-p,p' -P) + (™' -p')J^(p,<z) 

+{nY-Y')j^^{Y,q) - lAp.q) - lAp'.q) + L^q)] - {2k + ApY[{-mAp' +Y' - m') 

+p'Y^)V^{p,Y,q) + ip^+Y^ - 2m2)[/-(p,p',g) + 7-(q-p,p' -p) + {m^ - p^)r {p, q) 
+ {nY-Y')r{Y,q) - lAp,q) - I^P' ,q) + L'^iq)] + A[{-rY{p'' +Y^ - nY) 
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^ St 








-D12 




-Die ~ -DisV/'' y J: 


. r^':, = r^':,(y,g',p), 


^ -D18 








^ -D26 





+ (m2 - p2) JA.-(p, + (^2 _ p'2) ^) _ j^^^p^ _ j,^(^p'^ q) + ^^^-(5)]] } . (27) 

The rest of the diagrams can be obtained as follows: 

(28) 
(29) 
(30) 
(31) 
(32) 
(33) 

while the contributions F^*^, T^*^^, T^*^^, T'^^^ +^£1^2 vanish. Evaluation of the integrals completes the calculation 
of the 1-loop Compton scattering amplitude for off-shell external legs in arbitrary gauge and dimensions in SQED. 
Results for all the scalar and multi-indexed tensor integrals have been presented in the appendices. Through their 
most general decompositions in terms of the available vectors, the tensor integrals are converted into a series of scalar 
ones for their subsequent evaluation. Most of these results are new. We verify them to reduce to known expressions 
for simpler cases whenever possible. Moreover, as another confirmatory check, we verify the properties of the scalar 
pieces which reflect the symmetries of the original integrals. 

Once we know the full vertex at the one loop order, we can subtract from it the longitudinal part of the previous 
section to extract the unknown transverse coefficients Sn, -622, B23, -632 and -633 to 1-loop in perturbation theory. 
Although straightforward in principle, it is rather involved in practice. The outline is as follows. The 4-point Compton 
scattering vertex can be expanded in a natural basis as : 

r,,^=Co + Ci k^ky -t- C2 k^py + C3 p^ky + d k^p'^ + C5 p^fc^ + C& p^p'^ + C7 p^Pi. + Cs + Cg p,,p^,(34) 

where the coefficients Ci can be identified from the expressions (|19ll27p . Once the Ci are known, then comparing 
Eq. ([4]) and Eq. (j34|l we can write the unknown transverse pieces in Eq. ([4]) at the 1-loop level as follows : 

_ fc2 ^{k-p + p^ -P-P') ^ ^-P n ^ k^i-p'^+p-p' + k-p') 

k-p' „ k ■ p{-p'^ + p ■ p' + k-p') k-p'{k-p + p'^ ~P-P') ^ 

k-k' [k- k'Y ' ^ (fc • k'Y ' 

, k ■ p' {-p'"^ + p ■ p' + k-p') k ■ p{k ■ p + p^ - p ■ p') 

1 Cg + C7 + Cs + Cg 

i>22 
-B23 

-B32 

~ A{k-p'){~2k-p + k-p') + {k-p)^ + k'^{k^ + 2k-p-2k-p') ' ^ ' 

Thus the decomposition of the Compton scattering vertex is completely defined by Eqs. ([SllH]) and Eq. ([35|) . The 
longitudinal part has been determined non perturbatively, whereas the transverse part has been evaluated to 1-loop 
order. The perturbative transverse pieces can guide us to their possible non perturbative structures which should 
reduce to this expansion in the weak coupling regime. 



4(fc-p')(- 


-2k 


■ p + k ■ p') + {k -pY + k^{k'^ -\-2k - p- 


2k 


■P') ' 


1 




Cq ~ Cr + Cs — Cg 






4(fc.p')(- 


-2k 


■ p + k ■ p') + (k -pY -I- A:2(fc2 + 2k - p- 


2k 


■P') ' 


1 




2 C2 + 2 C4 — Cg + C7 + Cg — Cg 






4(fc-p')(- 


-2k 


■p + k ■ p') -\- {k ■ pY -\- k'^{k'^ -\-2k-p- 


2k 


■P') ' 


1 




2 C2 — 2 C4 + Cg + C7 — Cg — Cg 
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IV. CONCLUSIONS 



In this article, we provide a suitable basis to decompose the full Compton scattering vertex in SQED into its 
longitudinal and transverse (to one or both external photons) parts. We then employ the WFGTI which relates 3- 
point vertices to the 4-point ones in order to determine the longitudinal component of the Compton scattering vertex 
in an exact non perturbative fashion, written in terms of lower point functions. Recall that the WGTI-conserving 
model building for the three point interactions provides us with a reliable truncation scheme for the SDEs at the level 
of two-point functions. The work presented in this article provides an opportunity to lead us to the next order of 
approximation. One can now truncate the tower of SDE at the level of the 4-point function, again with a WFGTI 
conserving ansatz and consistently solve the coupled set of SDEs for the 2- and the 3-point functions in order to 
arrive at their more accurate solutions. A natural next step is to probe the non perturbative and vitally important 
structures of 4-point vertices in QCD through the generalized Slavnov- Taylor identities. 

As the name suggests, the transverse piece remains undetermined by the WFGTI. However, with the help of our 
complete 1-loop evaluation of this vertex, the unknown transverse pieces can be extracted to that order. It should 
serve as a guide for its possible non perturbative extensions. Any such attempt should reproduce this perturbative 
result in the weak coupling regime. 

There are further advantages of calculating Compton-scattering for massive electrons and off-shell external legs 
in arbitrary gauge and dimensions: (i) As the external legs are off-shell, one can study various cases of interest by 
putting the desired particles on shell. On the other hand, it can also serve as the internal part of more elaborate 
Feynman diagrams, (ii) Arbitrary gauge helps in checking the gauge invariance of the related physical observables 
and the gauge covariance properties of the Green functions, (iii) Substituting Z? = 4 — 2e and expanding in powers 
of e, we can study the 4-dimensional case. Interest has also been shown in lower dimensional SQED which can be 
simply projected out by an appropriate substitution of D. (iv) Explicit results of most of the integrals are new and 
can also be used in other calculations. These are the same integrals which are also encountered in spinor QED and 
QCD. Therefore, this study can provide a guideline for a similar calculation of the ghost-ghost-quark-quark, 4-gluon 
and 4-quark vertices in QCD, whose importance has been advocated in the introductory section. 



ACKNOWLEDGMENTS 



We are highly grateful to A. Davydychev for a very careful reading of our manuscript and suggesting to us several 
useful amendments, corrections and improvements. We acknowledge CONACyT, COECyT and CIC grants under 
projects 46614-1 and 94527, C8070218-4 and 4.10 respectively. 



APPENDIX I: Master Integrals 

In this appendix we summarize the results for various integrals that have appeared throughout the calculation. 
Although some of the integrals are basic and already known, we tabulate them all for the sake of completeness. All 
through the appendix, we use the simplifying convention Xq = (2/i7r^)A" for all scalar integrals with or without tilde. 



A. Scalar Integrals 



The 1- and 2-Point Integrals : 



We start from the list of all the scalar integrals up to 2-point integrals 



-z^^/2(m2)^/2-ir(^l 



2 ^ 

^2 



Kiq) = *7r^/2 (m2)^/2-2 r (^2 
Liq) = -z^^/2(m2)^/2-3r(^l 



(36) 



Kiq) = -*.^/^(m^)^/-^r(l - |) (2 - |, 1; |; ^ ) , (37) 



D\ p A, D D q' 



--)2Fi(3--,2;-;^) . (39) 
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Some of these scalar integrals can be compared with those reported elsewhere. 

• One can check that K{q) is the same as J(a, /?, 0, m) in Eq. (10) of reference [ll| for a = /3 = 1. 

• Similarly L{q) is the same as J{a, P,0,m) (as above) for a = 2 and (3 = 1. Moreover, using the recurrence 
relation (A. 17) of [l^. one can see that L{q) is not independent. It can be written in terms of T and K{q) as 
follows : 

^(9) = - (g2_^^2)2 [(D - 3) {q' + m^) K{q) - (D - 2)T] . (40) 

• Eq. (17) of [m in the case i^i = 1^2 = 1, P = q reproduces K{q). It can also be obtained from J(l, 1; mi, 771,2), 
Eqs. (A.l, A. 7) of [l^ for nii = m2 and x = y = rr? jq^ and after simplification. 



The 3-Point Integrals : 

We now come to the 3-point function of two variables fc and p : 



J(fc,p) 
J{k,p) 



.r(3-f) 



(^2)^/2-3^^__^^^ 



D 



D 

21-^1-^1-^1 2 
D. 
2 ' 



2 



3-f, 1,1,1 



fc2 (/c-p)2 p2 



771^ 



777^ 



777^ 



(41) 
(42) 



2x' 



-2(1) - 4){k~p)^ {ir? + (fcp)] /(fc,p) + 2(D - 3){k-p)^k{k-p) - {D - 2) 



{k-pf 



T 



(p2 - ra'^){k - p)2 - 2(p2 _ fc2^^2 r 

(p2 — 77i2)2 

(P - 777^)(fc - pf + 2(p^ - fc^)7n^ 
(fc2 — m2)2 



2(D - Z)p'K{p) - (D - 2) 



27772 



2(D-3)fc^if(fc) - (£'-2)- 



27(1,2 



(43) 



where x = "rn^Ot^ — P^)^ + ("^^ —k'^){'m? —p'^){k — pY. Note that the the integrals I{k,p) and J{k,p) are J2(l, 1, 1) and 
J2(l, 1, 2) of [6|. Generalized Lauricella function $2 arises for general powers of three propagtors z^i, 1^2 and t's (with 
the propagator corresponding to being massless, see [ill [igI . flTj and references therein, for definition, properties 
and symmetry relations) : 



$2 



^^123 ^ "J , i^l I 1^2 , ; "J ^1^1 
T'^23 



Xi,X2,X3 



Jj 1 JjH Jj 



2 -^3 



nd,j=0 



{1^123 - %]n + l + j) {vi;n + I) {v2]n+j) {1^3; I + j) - 
{§;n + l+j) iiy23;n + l + 2j) 



(44) 



where the Pochhammer symbol (z; 77) = T(z + n)/T(z) and h'i-^^i2,---,ij = I'ii + I'ia + ' ' ' + '^ij ■ Similarly, the generalized 
Lauricella function $3 ^ arises for general powers of three massive propagtors: i^i, 1^2 and 1^3, 



$3 



1^123 - T''^1''^2,J^3 
1^123 



Xi,X2,X3 



E 

n,l,j=0 



D . 



X'{ X2 X^3 ['^123 2 

77! n j! 



;n + l+j) {1^1,1+ j) {1^2; n + j) {1^3; n + I) 
ii^i23;n + l + 2j) 



.(45) 



We finally go on to the 4-point integrals. 



^ Note that the generahzed Lauricella function was first introduced in [29| . Some special cases, including those denoted as <1?2 and <1?3, 
were employed in lllll . Further examples can be found in [ll,[l3. 
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The 4-Point Integrals : 

Genuine 4-point integrals are eiven by Eq. (fT6|) and Eq. (fT7|) . We evaluate them using the Mellin-Barnes 
technique developed in [HI, [H, [13^ : 

U{j>,q,p') = j7r^/2(,„2)D/2-4£(i_li)£(i:_l) ^ 



^3 



4-— 1111— -1 

D o 
2 ' 



f_ (f_ ^ {q-vY {q - p'f {p-p'f 

919791 97 9 ' 9 



(46) 



3^3 



6r 



1,1,1,2,^-2 



D_ 

2 , 2 

^ 4 
2 



p2 ^2 ^,2 (g_p)2 (^_p,)2 (p_p,)2 



(47) 



The scalar integrals C/(p, g,p') and V(ji-, q,p'), given by Eqs. ITZ)) are genuine 4-point functions with three massive 
propagators, the most non-trivial integrals calculated in this paper. The function 3^3 is given by 



3^3 



1 1^1234 



^^1234 — -J, i^l, ^^2, t^3, ^^4; Y ~ ^^4 
-J, ^^123 



1 



X2 



n.l,r.w.k.j—0 



n\ U r! It;! kl j\ 



D 



;n -\- 1 -\- r ^ w -\- k ^ j) {ui^r -\- k j) {1^2]^^ -\- 1 -\- r) (^^3; n -\- w -\- j) (z/4; / + it; + A;) (-^ - 1^4; n + r + j) 



(48) 



{§;n + l + r + w + k + j) (1/123; 2n -(- Z -f 2r ?« + fc + 2j) 



This completes the calculation of the scalar integrals. 



B. Vector Integrals 

The more straightforward vector integrals which are functions of just one momentum variable are 



K^{q) 



2g2 



T-im' ~q')K{q) , 



-K{q) 



= ^ [iq'-m')L{q)+Kiq)] . 
Relatively more involved vector integrals are functions of two momentum variables which we take up now. 
The /^(fc,p) Integral : 
The most general description of I^{k,p) is : 



(49) 
(50) 
(51) 



Iti{k,p) = —[kf,lA{k,p)+Pf,lB{k,p)] 



(52) 



where 



lA{k,p) = -^|r - {k-p)]K^{p - k) + \p\k' - m') - k-p{p' - m')]h{k,p) +p'Koip) - {k ■ p)Ko{k) 
lB{k,p) = lA{p,k) , 



(53) 



^ Reference develops the method. Based upon it, the 4-point functions were considered (along with the arbitrary number of legs) 
in II6II for equal masses and in [TtII for arbitrary masses. 
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and = {k ■ pf - k'^p^. 
The I''ik,p) Integral : 
Similarly, we can write I^{k,p) as : 

. 2 

It.{k,p) = —[k^lA{k,p)+p^lBik,p)], (54) 

where lA{k,p) and Isik^p) are given by the same expressions as lA{k,p) and lB{k,p) in Eqs. (|53p with the ~ now 
appearing for each quantity, i.e., with the replacements Kq{p — k) Ko{p — k), Io{k,p) Io{k,p), Ko{p) — > Ko{p) 
and Ko{k) Ko{k). 

The J''{k,p) Integral : 

In the same way, the most general description of J^{k,p) is : 

- 9 

ITT 

Jfj^ik,?) = —[k,,JA{k,p)+Pf,JB(k,p)], (55) 

where JA{k,p) and JB{k,p) are given by the same expressions as lA{k,p) and lB{k,p) respectively in Eqs. (|53p with 
the replacements Ko{p — k) ^ Io{k,p), Io{k,p) Jo(fc,p), Kq{p) — > Lo{p) and A'o(fc) ^ Lo{k). 

We now treat the most complicated vector integrals which depend upon three independent momentum variables. 

The U^'{p,q,p') Integral : 

The U^{p,q,p') integral can be expressed only in terms of p^, and p'^. Therefore, 

■ 9 

Uf,{p,q,p') = —[pf,UA{p,q,p') + qt,UB{p,q,p')+p'f,Uc{p,q,p')], (56) 

where 

UAip,q,p) = 2[p'2(p . q)2 _ 2{p ■ q){p-p'){p' ■ q) + q\p ■ p'^ + p^p' ■ q)^ ' p'^q^)] ^ 
{[{p ■ P'W ■ q) - V'\V ■ q)]Io{p,p') + [p'\^ - ip' ■ qfW, q) 
+ [{p ■ q){p' ■ q) - q\p ■ p')]Io{p,q) + [p''\p ■ q) - p'\^ + q\p ■ p') 
-{p- q+p-p')p' -q+ip' ■ qf]Ia{q-p,p' - p) + [p''^q^{-p^ + P ■ q + P ■ p') 

-p ■ q{p'^{p ■q) + q^{p- p')) + p^{p ■ qf + m^{-p'^{p ■ q) +p'^q'^ 

-q\p ■p')+p -qip-q+p- p') - {p' ■ qf)]Uo{p, q,p')} , 
UB{q,P,p') = Ucip',q,p)^UA{p,q,p') . (57) 

The V''{p,q,p') Integral : 

The integral V^''(p, q,p') can also be expanded out as follows : 

. 9 
ZTT 

- —[pMp^q^p') + q^VB{p,q,p') + p'^Vc{p,q,p')] , (58) 

where Va{p, q,p') = VB{q,P-,p') = Vcip', q,p) are the same as C/a(p, q,p') = UB{q,P-,p') = Uc{p' , q,p) in Eqs. ([57]) with 
the following replacements : Io(p,p') Jo{p,p'), Io{p',q) Jo{p',q), Io{p,q) ^ Jo{p,q)- Io{q-p,p'-p) ^ Uo{p,q,p') 
and Uo(p,q,p') Vo{p,q,p'). 
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C. Tensor Integrals with Two Indices 



These integrals can also be categorized as regards the number of independent momentum variables they depend upon. 
For only one momentum dependence, we have 



1 



4:{D - 1 

1 



2{D - 2)T + {q^D - Am^) K{q) 



q^ d'^ 



Solutions for the integrals with growing complexity have been listed below separately : 
The I^"'{k,p) Integral : 

The tensor integral J^'' can be expressed in terms of the scalar integrals Kq, Ic, Id and Ie as follows : 



(59) 
(60) 



7-(2) ^ 1^ 

V 2 



^Ko{p - fc) + ( k^K ~ 9t^u-j^ ]lc + [ Pf^ku + kfj,p^ - 



D 



2(fc-p) 
' D 



Id + { Pt^Pv - 9tiuj^ ]Ie 



The coefficients /c Id and I^; in the above expressions are : 
D 



Ic{k,p) 



1 1 
2^ D 
2 1 
2 



[p — m )k ■ p 



1 1 

2^ D 



1 1 

2^ 2D 

Koip~k) 



{k^ - m^)p^ 



1 

2D' 



lA-—{p^-rn^)p'lB 



(k-p) 



lD{k,p) 



D 



2{D - 2)A2 



(fc^ - TO^)fc -p - [I ]ip- m^)k 



D 

/£ _ 1 

2 4 



2 



1 

D 

.2 ™2 



(61) 



(p — m )k ■ p 



1 _ ^ ) _ ^2)^2 



i^o(p-fc) 



(62) 



Note that the momentum dependence of Ia and Ib has not been displayed explicitly. It is lA^k^p) and lB{k,p). 
Moreover, Ic{k,p) and Isik^p) arc symmetric functions under the exchange of k and p. Therefore, 



lE{k,p) = IciP.k) . 

The Jf"'{k,p) Integral : 

The tensor integral J^" has the following decomposition in terms of the scalar integrals Io,Jc, Jd and Je 



(63) 



7(2) _ 1^ 



+ ( k^k^ - g^j,u'— ]Jc + [ P^ki, + k^p^, - g^^^ 



D 



D 



2 {k-p) 
D 



The coefHcients Jc, Jd and Je in the above expressions are : 
D 



Jc{k,p) 



{D - 2)A2 



1 1 

2^ 15 



{p — m )k ■ p 



1 1 

2^ 2D 



(fc2 - m^)p^ 



Jd + ( Pf^Piy - 9tiiyj^ ) Je 



J A ~ ^(P^ ^ ■n?)p^JB 



D 



1 



D 



2 



'4-^lB + ^Io-{l-i,]p.LUk)}, 



2D 



D 



D 



(64) 



Joik.p) 



D 



2(D-2)A2 



(fc^ - m^)fc -p - [I - — \{p^ - m^)k^ 



_(l_^)(fc2_,^2)p2 



2 



2 



(p^ — m?')k ■ p 



k.p-[i-^]e 



Ia 
2 



k-p-[i-%]p' 



2 



D' 



D 



-ik-p)Io +1-7^ [k.Kik) +P.L^oip)] 



(65) 
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The momentum dependence of Ia,Ib, Ja, Jb is lA{k,p), lB{k,p), JA(k,p), JB{k,p). Moreover, as is evident. Jc{k,p) 
and JE{k,p) are symmetric under the exchange of k and p, i.e., 



JE{k,p) = Jc{p, k) . 

The U''''{p,q,p') Integral : 

The general expansion of U^" in terms of the scalar integrals Iq,Ud, Ue ,Uf, Ug, Uh and Ui is as follows : 



Uf,i,ip,q,p') = — 



where 



Ud{p, q,p') 



^lo{q -p,p' -p)+ ( Pf_,p,y - g^u^ Wd+[ Pp^qu + q^Pu - g^J,u 



D 



+ ( Pt^Pu + Pp.Pi' - 9i^y — — 
2(y-<z) 



D 



Up + ( qtiqu - gt^iy-^ ]Ug 



2 (P ■ q) 
D 



Ue 



+ 1 qtJ.p'u + p'^iiu - fJiiu 



D 



Uh 



Pt^Pu -9t^i^jj ]Ui 



2{D - 3) [p^ip' ■ qf - 2{p ■ p'){p ■ q){p' ■q) + {p- p'fq^ + p'"^ {{p ■ qY - p'^q^)] 
{ [{D - 2)p\p' ■ qf -{D- i){p'^p -q + p- p'q^)p' ■ q + p'^ {{D - 3){p ■ p' + p ■ q) 
-{D - 2)p^)q^ + m^i-{D - 2){p' ■ qf + (D - 3)(p ■ p' + p ■ q)p' ■ q 
-{D- 3)p'^p -q+iD- 2)p'^q^ -{D- 3)p ■ p'q^)] Ua{p, q,p') 
+ [im'-q'-){p''q^~ip'-qf)]UB{p,q,p') 
+ [{p''-m'){{p'-qr^p'\')]Ucip,q,p') 

+ [{D - A){p' ■ qf -{D- i){p-p'+p- q)p' -q+iD- i)p'^p ■ q 
-{D- A)p'\^ + {D- 3)p • p'q^] To{q ~p,p'~ p) 

-iD-3) [{p' -qf-ip-p'+p- q)p' ■q + p'Hp-q-q^)+p- p'q^] X 
[lAiq-p,p' -p) + lB{q-p,p' -p)]+iD-3) [{p' ■q){p-p')- p'^p ■ q] X 
I A [p, p') + {D-i)[{p' ■q){p-q)-p- p'q^] I a {p, q) } , 



(66) 



(67) 



(68) 



Uh{p, q,p') 



2{D - 3) [p2(p' . g)2 _ 2{p-p'){p ■ q){p' ■ q) + {p-p'fq^ + p'^ ((p • qf - p^q^)] 

{ - [(771^ -p^){{p' ■q)p'^ - {p-p'){p-q))]UA{p,q,p') 

+ [{D - 3)p ■ q{p'^p -q-p'- qp") + (-p'((i? - 3)p'2 +p'-q + ?,p-p' 

—D{p' ■ q + p ■ p )) - {D - \)p ■ pp ■ q)q^ + m^(p' • q{-Dp^ + p^ 

+ {D -i)p.q)-p. q{-Dp -p' +p-p' + (D-3)p-q) + {D-3){p^ -p- p')q'-)] 

~p ■ p'iiD - 3)p ■ p' - Dp ■ q + p ■ q)) + p'^iD - l)p' ■ qp" 

+ {D -3){p-q- q^)p^ -{D- l)p ■ p'p ■q) + {D- 3)p ■ p' {p ■ p' q^ - p' ■ qp^)] 
Uc{P,q,P') , r. ,w ^ , 217/ ' \ 

\- [{p ■ p ){p ■ q) - p -qp \ io{q-p,p -p) 

\ 2 I I , / 2 2 2 I / 21 

n ^ — [p- q -p-qp-q-p-pp-q + p-qp -p q +p-pq \ 

~ D — 3 

lA{q -p,p' -p)-\ — [p' ■ q[p^ - p-p')+p- p'[p -p' -p-q) + 

p" {p-q-p^}] Ib {q -p,p'-p) + [p'V ~{P-P'f] Ib {P, P') 

+ ^2'^ [P^g^ - {P- qf]lB{p, q) + ^ ^ ^ [p' ■ qp - q-p-p'q^]lA{p' ,q) 



2 

D-3 



D — 3 

H — 2 — 'IP'P' ^ P'^P ■ ^b{p', q)] ■ 
UG{q,P,p') = Ui{p',q,p) ^ UD{p,q,p') , UE{p,q,p') ^ Uh{p ,q,p) UF{p,q,p) = UH{q,P,p) 



The V>'''{p,q,p') Integral : 

Similarly to the previous case, the V^^'^ integral can be expanded in terms of Uq.Vd: Ve ,Vf, Vg, Vh and Vi as 



V^l\p,q.p') = ^ 



^Uaip, q,p') + (^Pt^P,, - gf^i'j^^ Vd + (^Pt^qi^ + qiiPu - .9^1.^^^^ I Ve 
+ {pt^p'u + p't_iPv - g^-^ ^^^j)^ + (^9^9-^ - Spi'lj^ 



where the independent scalar integrals are : 

VD{p,q,p') = 



2{D - 3) [p2(p' . qf - 2{p ■ p'){p ■ q) ip' ■q) + {p- p'fq'^ + p'^ {{p ■ qf - p^q^)] 
{ [{D - 2)p\p' ■ qf -{D- 3)(p'2p .q+p. p'q^)p' ■ q + p'^{{D - 3){p ■ p' + p ■ q) 
-{D - 2)p2)q2 +rn^-{D- 2)(p' • qf + {D ~ 3){p ■ p' + p ■ q)p' ■ q 
-{D - 3)p'^p -q+iD- 2)p'\^ -{D- 3)p ■ p'q^)] Va{p, q,p') 
+ [(m2 - q') {P'V ip' ■ qf)]VBip,qy) + [{p" - m^) ((p' • qf - p'V) ] 
Vcip, q, p') + [{D - 2){p' ■qf-{D-3)ip-p'+p- q)p' -q + iD- 3)p'^p ■ q 
-{D- 2)p'\^ + {D^ 3)p ■ p'q^] Ua{p. q,p') + [{p' ■ qf - p'\^] [Ub{p. q,p') 
+Ucip, q,p')] + (£> - 3) [p -qp-p - p'^p ■ q\ Ja{p,p') 

+ {D-S)[p' -qp-q-p- p'q^] Ja{p, q) + 2 [p'\^ - {p' ■ qf] Ua{p, q,p')} , 



VH{p,q,p') 



2{D - 3) [pHp' ■ qf - 2{p-p'){p ■ q){p' ■ q) + {p.p'fq^+p'^ ((p • qf-p^q^)] 
{ - [(m^ - p2) (p' -qp^ -p- p'p ■ q)] Va{p, q, p') 

+ [{D - 3)p • q{p'^p -q-p'- qp") + ( - 3)p'2 +p' ■ q + 3p ■ p' 

-D{p' -q + p- p')) -{D- l)p • p'p • q)q^ + m'^{p' ■ q{-Dp^ + p^ + {D - 3)p • q) 

p-q{-Dp-p'+p-p' + {D- 3)p ■q) + {D-3){p'-p- p')q^)] ^^^^^ ^'^ 



[m'^ip' ■ q{-Dp^ +p^ + {D- 3)p • p') + {D- 3)p'\p'^ ~ p ■ q) - p ■ p'{{D - 3)p ■ p' 
-Dp-q+p-q))+ p'\{D - l)p' -qp" + {D-i){p-q- q^)p^ - {D - l)p • p'p • q) 



+ {D-i)p-p'{p-p'q^ -p' -qp")] 



Vc{p, q,p') 



+ [p' -qp^ - p- p'p ■ q] Ua{p, q, p') 



[{D - l)p' • qp^ - {D — 3)p' ■ qp- q+p ■ q{-Dp ■ p + p ■ p + [D - 3)p • q) 



+ (D - 3){p ■ p' - p^)q 



21 UB{p,q,p') 



+ [{D ~ l)p' ■ qp^ - (D - 3)p' ■ qp ■ p' 



+ {D - 3)p'2(p .q--p^)+p. p'{{D -3)p-p' ~Dp-q + p-q)] 



Uc{p, q,p') 



{D 3) r ,2 2 / /\21 T / l\ , 3) r 2 2 f \21 



■[p'V-{p-p'f]Mp,p) + 



■[p\^ ~iP-qf]Mp,q) 



+ [p' .qp.q-p. p'q^] Ja{p' ,q) + ^ [p' ■ qp ■ p' ~ p'^p ■ q] Jb{p\ q) 

+2[p ■ p' p ■ q - p' ■ qp^]Uo{p,q,p')] , 
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VG{q,p,p') ^ Vi{p\q,p)^VD{p,q,p') Veip, q,p') ^ Vnip' , q,p) Vf{p, q, p') ^ Vniq, p,p') ■ (74) 
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